Mon. Not. R. Astron. Soc. 000, [L>?? (2012) 



Printed 12 February 2013 



(MN LMfeX style file v2.2) 



Proper motions and membership probabilities of stars in 
the region of open cluster NGC 3766/ 



R. K. S. Yadav 1 ! Devesh P. Sariya 1 ' 2 , R. Sagar 1 



1 Aryabhatta Research Institute of Observational Sciences (ARIES), Manora Peak, Nainital - 263129, India 
^ 1 2 School of Studies in Physics & Astrophysics, Pt. Ravishankar Shukla University, Raipur-492 010 (CG), India 

Accepted Received ; in original form 



o 

6 



> 
oo 
in 
in 

(N 

o 



is 



ABSTRACT 

Relative proper motions and cluster membership probabilities (P^) have been derived 
for ~ 2500 stars in the field of the open star cluster NGC 3766. The cluster has 
been observed in B and V broadband filters at two epochs separated by ~ 6 years 
using a wide-field imager mounted on the WFI@ES02.2m telescope. All CCD frames 
were reduced using the astrometric techniques described in Anderson et al. (2006). 
The proper motion r.m.s. error for stars brighter than V ~ 15 mag is 2.0 mas yr _1 
but it gradually increases up to ~4 mas yr _1 at V ~ 20 mag. Using proper motion 
data, membership probabilities have been derived for the stars in the region of the 
cluster. They indicate that three Be and one Ap stars are member of the cluster. The 
reddening E(B — V) = 0.22 ±0.05 mag, a distance 2.5±0.5 kpc and an age of ~ 20 Myr 
are derived using stars of P u > 70%. Mass function slope x = 1.60 ± 0.10 is derived 
for the cluster and cluster was found to be dynamically relaxed. Finally, we provide 
positions, calibrated B and V magnitudes, relative proper motions and membership 
probabilities for the stars in the field of NGC 3766. We have produced a catalog that 
is electronically available to the astronomical community. 
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1 INTRODUCTION 

The similar age of the stars in open clusters make them vital 
objects for the study of stellar evolution. However, for such 
studies, it is essential to separate the field stars from the 
cluster members. Proper motions (PMs) of the stars in the 
region of open clusters provide a unique possibility of get- 
ting membership information of the stars. Proper motions 
in combination with radial velocities also provide valuable 
information on the kinematic parameters of the Galactic ro- 
tation curve and the distance to the Galactic center. They 
are also useful to study the formation and evolution of the 
open cluster system. 

The open cluster NGC 3766 (q 2 ooo = ll h 36 m 14. s 0; 
feooo = -61°36'30"; I = 294? 12; b = -0?03) is located in the 
Carina complex of the Milky Way. The cluster has been the 
target of many photometric studies (Ahmed (1962); Yilmaz 
(1976); Shobbrook (1985), (1987); Moitinho et al. (1997); 
Tadross (2001); McSwain & Gies (2005)) with the aim to de- 



* Based on observations with the MPG/ESO 2.2m, located at 
La Silla and Paranal Observatory, Chile, under DDT programs 
164.O-0561(E) and 077.C-0188(B). 
f E-mail: rkant; devesh; sagar@aries.res.in 



termine the various physical parameters of the cluster. The 
most remarkable attribute of NGC 3766 is the presence of 
largest number of Be stars yet known in any Galactic clus- 
ter. Among the 37 brightest stars, 12 show emission lines 
(Sher (1965); Schild (1970); Lloyd Evans (1980)). Funda- 
mental parameters of B stars in NGC 3766 was studied by 
Aidelman et al. (2012) using spectroscopic data. 

A proper motion study of NGC 3766 was made by 
Baumgardt et al. (2000) using eight stars taken from Hippar- 
cos catalogue and derived the mean absolute proper motion 
/j. a cos(5) = —7.21 ±0.34 masyr -1 and /is — 1.33 ±0.35 mas 

yr" 1 . 

In spite of extensive photometric studies of this cluster, 
there is a lack of study providing proper motions and mem- 
bership probabilities (P M ) of the stars in the wide-field re- 
gion of open cluster NGC 3766. With the photometric data, 
membership information is desirable to minimize field star 
contamination. The archival wide-field multi-epoch CCD ob- 
servations taken with the WFI@2.2m telescope provide an 
opportunity to derive precise PMs with a time-span of only 
a few years, and is deeper by several magnitudes than pre- 
vious photographic surveys (Anderson et al. 2006; Yadav et 
al. 2008; Bellini et al. 2009; Sariya, Yadav & Bellini 2012). 

The main purpose of the present study is to estimate the 
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Table 1. Description of the WFI@2.2m data sets used in this 
study. First epoch data was observed on Feb 27, 2000 while second 
epoch data was observed on April 9, 2006. 



Filters Exposure Time Seeing Airmass Limiting Mag. 
(in seconds) 

(First epoch) 

B 2x30; 2x240 1".5 1.25 ~21 

V/89 1x30; 2x240 1".5 1.33 ~20 

(Second epoch) 

V/89 4x50 1".2 1.35 ~20 



fundamental parameters of the cluster using stars chosen ac- 
cording to their derived membership probabilities {P^)- The 
of different sources in NGC 3766 is also discussed. We 
also derived luminosity function (LF), mass function (MF) 
and mass segregation of the cluster using stars of P M > 70%. 
The PMs, P^ and photometric B and V magnitudes are 
provided to the astronomical community for the follow-up 
studies. 

Observational data alongwith reduction procedures and 
determination of proper motions are described in the next 
Section. Colour-magnitude diagrams using proper motions 
are described in Sec. [3] In Sect.[4]we present cluster member- 
ship analysis. Sect.[5]is devoted to derivation of fundamental 
parameters of the cluster. Luminosity and mass function are 
described in Sec.[6]while SecQis devoted to the mass segre- 
gation study. In Sect.[H]we describe our catalogue and finally, 
in Sect. [9] we present the conclusions of present study. 



2 OBSERVATIONAL DATA AND 
REDUCTIONS 

CCD data for NGC 3766 were collected with the wide-field 
imager camera (WFI) mounted on 2.2m ESO/MPI tele- 
scope at La Silla Chile. Using these archival dat£0 of two 
epochs, proper motions were computed. The first epoch con- 
sists of four images in B filter and three images in V filter 
taken on 27 th Feb 2000, while second epoch have four im- 
ages in V filter taken on 9 th April 2006. The observational 
log is listed in Table Q] 

The WFI@2.2m consists of eight 2048x4096 EEV 
CCDs with 0'.'238 pixel -1 resulting in a total field-of-view 
34' x 33'. Images used in the analysis were taken between 
1".2 — 1".5 seeing condition and between 1.25-1.35 airmass. 
Long and short exposures were acquired to map the brighter 
as well as fainter stars of the cluster. 

2.1 Astromeric and Photometric Reductions 

For the reduction of mosaic CCD images, we adopted the 
procedure described in Anderson et al. (2006, Paper I). The 
procedure include de-biasing, flat-fielding and correction for 
cosmic rays. To get the position and flux of the stars in the 
image, we made an array of empirical Point Spread Func- 
tions (PSFs) for each image. An array of 15 PSFs in each 
2048x4096 pixels chip (3 across and 5 high) as shown in Fig. 

1 http: //archive. eso.org/eso/eso_archive_main. html. 
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Figure 1. Plot of rms of the residuals around the mean magni- 
tudes in B and V and as a function of V magnitudes. 

3 of paper I was considered because PSFs changes signifi- 
cantly with position on the detector. In total, 120 PSFs have 
been considered for entire field of view (8192x8192 pixels). 
These PSFs are saved in a look-up table on a very fine grid. 
To select suitable stars for the PSFs, an automatic code was 
developed (see Paper I). An iterative process is designed to 
work from the brightest down to the faintest stars and find 
their precise position and instrumental flux for B and V 
exposures. 

In paper I, it is shown that WFI@2.2m has a large 
geometric distortion, i.e. pixel scale is changing across the 
field of view (see Paper I). To derive the correction for ge- 
ometric distortion,we parametrized the distortion solution 
by a look-up table of corrections for each chip that cov- 
ered each 2048x4096 pixel chip, sampling every 256 pixels. 
This resulted in a 9x 17 element array of corrections for each 
chip. At any given location on the detector, a bi- linear in- 
terpolation between the four closest grid points of the look- 
up table provided the corrections for the target point. The 
derived look-up table may have a lower accuracy on the 
edges of a field because of the way the self-calibration frames 
were dithered (see Paper I). An additional source of uncer- 
tainty is related to a possible instability of distortions for 
the WFI@2.2m reported earlier. This prompted us to use 
the local-transformation method to derive PMs. Detailed de- 
scriptions about the distortion solution is given in Paper I. 

In the local transformation approach a small set of local 
reference stars is selected around each target object. It is 
advantageous to use pre-selected cluster members to form 
a local reference frame because of a much lower intrinsic 
velocity dispersion among the cluster members. Then, six- 
parameter linear transformations are used to transform the 
coordinates from one frame into another, taken at different 
epochs. 

A general linear transformation from one frame (x,y) 
into another (u, v) has 6 parameters, and can be put into 
the following form: 
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The parameters A, B, C and D are the linear terms, and 
xo, yo, uo and vq are the constant terms. The equation look 
like there are 8 free parameters in the equations, but we can 



Astrometry of NGC 3766 3 



choose the zero-point in one frame arbitrarily, so only one of 
the two offsets is actually solved for. The easiest way to solve 
for the offsets is to adopt the centroid of the star list in the 
first frame as the offset for that frame, (xo,yo). Then, the 
centroid in the other frame will be the least-squares solution 
for the other offset, (uq,vo). 

The residuals of this transformation are characterizing 
relative PMs convolved with measurement errors. In essence, 
this is a classical "plate pair" method but extended to all 
possible combinations of the first and second epoch frames. 
Relative PM of a target object is the average of all dis- 
placement measurements in its local reference frame. The 
last step is to estimate the measurement errors from intra- 
epoch observations where PMs have a zero contribution. A 
complete description of all steps leading to PMs is given in 
Paper I. 

2. 1.1 Photometric Calibration 

In order to transform instrumental B and V magnitudes 
into the standard Johnson-Cousin system, we used the data 
published by Moitinho et al. (1997). Total 92 stars were used 
to calibrate B and V magnitudes. 

We derived photometric zero-points and colour terms 
using the following transformation equations: 
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Figure 2. Plots of proper motions and their rms error in 
mas yr versus visual magnitude. 



-Bstd = -Bins + Cb * (Bins — Vins) + Zb 
V st d = Vi ns + C v * (B ins — Vins) + Z v 

where the subscript "ins" means instrumental magni- 
tudes and "std" stands for secondary standard magnitudes. 
Cb and C v are the colour terms while Zb and Z v are the 
global zero-points. The quadratic colour terms are negligi- 
ble. The values of colour terms are 0.29 and 0.16, while the 
zero-points are 24.39 and 23.79 for B and V filters respec- 
tively. 

In Fig. [T] we show the rms of the residuals around the 
mean magnitude for each filter as a function of the V magni- 
tude. The photometric standard deviations have been com- 
puted from multiple observations, all reduced to the com- 
mon photometric reference frame in the chosen bandpass. 
The stars brighter than ~ 12.0 mag in V have higher dis- 
persion because of saturation. On average, photometric rms 
are better than ~0.02 mag for stars between 12-18 mag and 
~ 0.04 mag between 18-20 mag in V filter. 

2.1.2 Astrometric Calibration 

The X, Y raw coordinates of each star in each frame were 
corrected for geometric distortion using the look-up table 
provided in Paper I, brought into a common reference frames 
by means of six-parameter linear transformations, and aver- 
aged. To transform the averaged X and Y coordinate into 
RA and Dec of J2000, we used the online digitized sky ESO 
catalogue in skycat software, as an absolute astrometric ref- 
erence frame. Thanks to our accurate geometric-distortion 
solution and a reasonable stability of the intra-chip posi- 
tions, that it was possible to apply a single plate model 
involving linear and quadratic terms and a small but signif- 
icant cubic term in each coordinate. This solution also ab- 
sorbs effects caused by differential refraction. The standard 
error of equatorial solution is ~ 30 mas in each coordinate. 



2.2 Determination of proper motions 

Proper motions were computed using ^-filter images to min- 
imize colour-dependent terms in our analysis. Moreover, our 
geometric-distortion solution provides the smallest residuals 
with V images (Paper I). A total of three images for first 
epoch and four images for the second epoch were used in 
proper motion determination. 

First, we selected a sample of probable cluster members 
using the V vs. (B — V) CMD. Selected stars are located on 
the main-sequence in the magnitude range 8.0< V < 18.0 
mag. These stars define a local reference frame to transform 
the positions of a given first-epoch image into positions of 
a second-epoch image. By adopting only stars on cluster 
sequence having proper motion errors <3.0 mas yr" 1 , we 
assure that PMs are measured relative to the bulk motion of 
the cluster. To minimize the effect of uncorrected distortion 
residuals we used the local transformation approach based 
on the closest 25 reference stars on the same CCD chip. No 
systematics larger than random errors are visible close to 
the corners or edges of chips. 

We iteratively removed some stars from the preliminary 
photometric member list that had proper motions clearly 
inconsistent with cluster membership, even though their 
colours placed them near the fiducial cluster sequence. The 
distribution of proper motions rms error is presented in Fig. 
[2] The precision of the proper motion measurement is bet- 
ter than 3 mas yr -1 upto 14 mag in V. They are gradually 
increasing upto 4 mas yr -1 for V = 20 mag. 

2.2.1 Comparison with UC AC 3 proper motion data 

UCAC3 catalogue (Zacharias et al. (2010)) provides the ab- 
solute proper motions of the stars in the region of NGC 
3766. There are 94 common stars found to be brighter 
than V — 19.0 mag. To compare our proper motions with 



4 Yadav et al. 



UCAC3 data 
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Figure 3. Top panels: Vector-point diagrams of common stars 
relative to the cluster mean motion for UCAC3 (left) and our 
catalogue (right). (Bottom panels:) p. a CosS and ps as a function 
of X and Y in WFI pixels. Gaps in X and Y axis are the inter- 
chip space present in the mosaic CCD. 



UCAC3, we changed the UCAC3 proper motions to rela- 
tive proper motions. For this, we subtracted the mean abso- 
lute proper motion of the cluster (pi a cos8 — —3.4 mas yr -1 , 
Us — —0.9 mas yr" 1 ) from individual proper motion of stars. 
Fig. [3] shows the comparison of our PMs with those of the 
UCAC3 catalogue. Top-left panel shows the vector-point di- 
agrams (VPDs) of UCAC3 stars while top-right panel shows 
the VPD of our measurements. A concentration of stars 
around (0, 0) mas yr _1 is seen in both VPDs. Our proper 
motions distribution is tighter than the UCAC3 distribu- 
tion. This is because our data is more precise than UCAC3 
data. 

In the lower panels of Fig. [3] we show the proper motions 
in Hc,cos5 and pis as a function of X and Y coordinates. 
There is no clear systematic trend seen in this diagram. 



3 CLUSTER CMD DECONTAMINATION 

PMs pt a cos5 and (is are plotted as VPD in the top pan- 
els of Fig. 3] Bottom panels show the corresponding V vs. 
(B — V) CMD. Left panels show all stars while middle and 
right panels show the probable cluster members and field 
stars. A circle of 4 mas yr _1 around the cluster centroid in 
VPD of proper motions defines our membership criterion. 
The chosen radius is a compromise between losing cluster 
members with poor proper motions and including field stars 
sharing the cluster mean proper motion. The shape of clus- 
ter members PM dispersion is circular, providing that our 
PM measurements are not affected by any systematics. The 
CMD of most probable cluster members is shown in the 
lower middle panel. Main-sequence of the cluster is clearly 



separate out, which demonstrates the power of PMs derived 
in this study. The lower right panel represents the CMD for 
field stars. Few cluster members are also visible in this CMD 
because of their poorly determined proper motions. 

Fig. [5] shows (B — V), V CMD with VPD in various 
magnitude bins. We divided the CMD into five magnitude 
bins and in each bin, we adopted different selection crite- 
ria to identify cluster members, which were more stringent 
for stars with more reliable measurements from data of high 
signal-to- noise ratio, and less restrictive for stars with less 
precise measurements. These stars have a proper-motion er- 
ror of < 3.0 mas yr _1 . For each magnitude bin, we consid- 
ered as cluster members those stars with a proper motion 
within the circle shown in the middle column of Fig. [5] The 
right side of Fig. [S] shows the color-magnitude diagram for 
stars assumed to be cluster members. This figure shows that 
the separation of brighter cluster members from the field 
stars is clearly visible while fainter members are not clearly 
separated out. The reason may be that proper motions for 
fainter stars are not determined accurately. The available 
archive images are again not sufficiently deep to derive re- 
liable proper motions of fainter stars. The proper motions 
derived in this paper are also not sufficiently accurate to 
study the internal motion of the cluster. 



4 MEMBERSHIP PROBABILITIES 

The VPD for the proper motions derived in this study is 
shown in Fig|3] Present proper motion leads to the compact 
appearance of the cluster stars in VPD. Membership deter- 
mination based on proper motions is useful for further as- 
trophysical studies in the region of cluster. The fundamental 
mathematical model set up by Vasilevskis et al. (1958) and 
the technique based upon the maximum likelihood principle 
developed by Sanders (1971) for membership determination 
have been continuously used and refined. 

An improved method for membership determination of 
stars in clusters based on proper motions with different ob- 
served precisions was developed by Stetson (1980) and Zhao 
& He (1990). Zhao & Shao (1994) then added the corre- 
lation coefficient of the field star distribution to the set of 
parameters describing their distribution on the sky. 

A tight clump at jj, a cos8— pis =0.0 mas yr _1 in Fig. [4] 
is representing the cluster stars and a broad distribution 
is seen for the field stars. Fig. [5] indicates that the field 
stars centroid for panel 3 (middle column) is at a different 
location than for the last panel. Therefore, we consider two 
groups (V <= 15 mag and V > 15 mag) for membership 
determination. To determine the membership probability, 
we adopted the method described in Balaguer-Nunez et al. 
(1998). This method has already been used for uj Centauri 
(Bellini et al. (2009)) and by us for NGC 6809 (Sariya et 
al. (2012)). In this method first we construct the frequency 
distribution of cluster stars (<f>c) an d field stars (</>/)■ The 
frequency function for the i star of a cluster can be 
written as follows: 



4>" = 



i r (fx 



2 



+ 



where pi xi and pi y i are the proper motions of the i star 
while fi xc and fj, yc are the cluster's proper motion center. a c 
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Figure 4. (Top panels) Proper motion vector-point diagram. Zero point in VPD is the mean motion of cluster stars. (Bottom panels) 
Calibrated V vs. (B — V) CMD. (Left) The entire sample; (Center) stars in VPD with proper motions within 4 mas yr _1 around the 
cluster mean, (right) Probable background/foreground field stars in the direction of NGC 3766. All plots show only stars with proper 
motion a smaller than 3.0 mas yr — 1 in each coordinate. 



is the intrinsic proper motion dispersion of cluster member 
stars and (e x i,e y i) are the observed errors in the proper 
motion components of i th star. The frequency distribution 
for i th field star is as follows: 



while /j, x f and n y f are the field proper motion center. e xi and 
t v i are the observed errors in proper motions component and 
a x f and a y f are the field intrinsic proper motion dispersions 
and 7 is the correlation coefficient. It is calculated by 



0/ = 



2,V(l- 7 2)^ (CT 2^ +£ 2 



27<Mxi-Mx/)(Ml/i-M; / /) , (V-yi-V-yfY 



Mo 



Vf 



2 ) 



exp{- 
]} 
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where and jj, y i are the proper motions of i star 



7 = 



We did not consider the spatial distribution of the 
stars due to small observed field. To define the distribution 
function <ff c and <fff we used stars having proper motion 
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Figure 5. (Left:) CMD for stars with proper motions. (Middle:) VPD for the same stars in corresponding magnitude intervals. A circle 
in each plot shows the adopted membership criterion. (Right:) CMD for stars assumed to be cluster members. All plots show only stars 
with proper motion cr smaller than 3 mas yr — 1 in each coordinate. 



error better than 3 mas yr -1 . As expected, in the VPD, 
the center of cluster stars is found to be at x c = 0.0 
mas yr -1 and y c = 0.0 mas yr -1 . Our proper motion 
data set could not determine intrinsic proper motion 
dispersion (<j c ) for cluster stars. Assuming a distance of 
2.5 kpc (present estimate) and radial velocity dispersion 
1 kms' 1 for open clusters (Girard et al. (1989)), the 
expected dispersion in proper motions would be ~ 0.08 
mas yr -1 . Therefore, we adopted a c = 0.08 mas yr" 1 . For 
field stars V <= 15 mag, we have: xj = —4.9 mas yr _1 , 
yf — 2.7 mas yr~ , a x f = 5.7 mas yr~ and a y f = 4.0 



mas yr 1 , while for V > 15 mag: Xf = —5.4 mas yr 1 , yj 
= 0.5 mas yr -1 a x f = 6.6 mas yr -1 and a y f — 5.3 mas yr -1 . 

The distribution of all the stars can be calculated as 

<f>=(n c ■ <j>l) + (n f . (f>" f ) 

where n c and n/ are the normalized number of stars 
for cluster and field (n c + n/ = 1). For the first group 
(V <— 15 mag), n c = 0.64 and n/ = 0.36 while for the 
second group (V > 15 mag), n c = 0.25 and n/ = 0.75. 
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Figure 6. Membership Probability P M as function of the V mag- 
nitude, for all the stars in our catalogue. At V ~ 18 mag and 
fainter, P M diminishes as a result of increasing errors in the PMs. 
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Therefore, the membership probability for i 4 star is 
P*(i) = W 

A good indicator of cluster and field separation is the 
membership probability. It is plotted as a function of mag- 
nitude in Fig. [6] As seen in this plot, high membership prob- 
ability (Pp > 90%) extend down to V ~ 17 mag. At fainter 
magnitudes the maximum probability gradually decreases. 

Fig. shows the histogram of membership probabilities 
for 2468 stars in the cluster region. It shows a clear sep- 
aration between cluster members and field stars. This also 
indicates that the proper motion membership determination 
made using the method mentioned above is effective for the 
cluster NGC 3766. We find 818 stars with membership prob- 
abilities higher than 70%. 



4.1 Effectiveness of membership determination 

During the observations of the cluster we can not avoid 
the contamination by background and foreground objects. 
The effectiveness of our membership determination can 
be estimated using membership probabilities as described 
in Shao & Zhao (1996). The effectiveness of membership 
determination is estimated as: 



Figure 7. Histograms of membership probabilities of all the stars 
in the region of the cluster NGC 3766 



E : 



Nj2 N =1 lP(i)xg-P(i))] 



where the bigger E is, the more effective the member- 
ship determination is. 

So we determine the effectiveness of membership deter- 
mination as 0.65 for NGC 3766. It is shown in Fig. 3 of Shao 
& Zhao (1996) that the effectiveness of membership deter- 
mination of 43 open clusters are from 0.20 to 0.90 with peak 
value 0.55. Compared with the work by Shao & Zhao (1996), 
we can see that the effectiveness of membership determina- 
tion for NGC 3766 is significantly high. 



5 2.5 




3 4 5 

Radius (arcmin) 



Figure 8. Radial surface density distribution of the stars with 
P M > 70% in NGC 3766 region. 



4.2 Membership probability of published Be and 
Ap stars 

Present catalogue can be used to assign membership prob- 
ability of various types of stars in this cluster. We found 
three Be stars and one Ap star in our catalogue according 
to the categorization in WEBDA. Membership probabilities 
of these stars are listed in Table [2] The value of P M >= 90% 
for the listed stars indicate that they are member of the 
cluster. 
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Table 2. Membership probabilities of Be and Ap stars in the 
region of NGC 3766. 
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5 FUNDAMENTAL PARAMETERS 

5.1 Surface density distribution 

Fig. [8] shows the radial density distribution of stars having 
membership probability P M > 70%. It is easily seen from 
the figure that the surface density distribution of the stars 
with radial distance (r) from the field center greater than ~ 
5 arcmin is quite flat, roughly being p ~ 2 stars per square 
arcminutes, which is due to field stars with the same proper 
motion as NGC 3766. So we can reasonably choose 5 ar- 
cmin as the cluster radius. The cluster radius 4'. 6 derived 
by Moitinho et al. (1997) with stars V < 17 mag is lower 
than the present estimate. If the distance of the cluster de- 
termined by us is adopted (D = 2500 pc), the corresponding 
linear radius is 3.6 pc. 

5.2 Interstellar reddening, distance and age 

In order to determine the reddening and distance of the 
cluster, the (B — V), V CMD was plotted in the left panel of 
Fig. [9] for the stars P M > 70%. The reddening and distance 
is estimated by fitting the ZAMS given by Schmidt-Kaler 
(1982) to the observational CMD. The best coincidence of 
ZAMS line with CMD of the cluster was achieved at the 
value of color excess E(B — V) = 0.22 ± 0.05 mag. Our de- 
rived value of reddening agrees fairly well with the value 
E(B - V) = 0.20 ± 0.10 mag estimated by Moitinho et al. 
(1997). The best fit ZAMS provides the value of apparent 
distance modulus (v — M v ) = 12.75 ± 0.05 mag, which cor- 
responds to a distance of 2.5±0.5 kpc. For this cluster, our 
estimated value of distance is similar to the value 2.2±0.3 
kpc derived by Moitinho et al. (1997). 

The age estimate of the cluster is based on visually fit- 
ting the theoretical isochrones described in Girardi et al. 
(2000). Isochrones of log(age) = 7.20, 7.30 and 7.40 are 
shown with continuous lines in the right panel of Fig. [9] 
The superposition of various isochrones on the CMD indi- 
cates that CMD is represented well with an isochrone of 
about 20 Myr. The present age estimate is not too differ- 
ent from the value 21 Myr derived in the previous studies 
(Moitinho et al. (1997), Mermilliod & Maeder (1986)). 



6 LUMINOSITY AND MASS FUNCTION 

In order to estimate the luminosity function (LF), histogram 
of the number of stars were constructed with 1.0 mag inter- 
vals in M v using the stars of P M > 70% and shown in Fig 
1101 The interval of 1.0 mag was chosen so that there would 
be enough stars per bin for good statistics, as well as to 
provide a reasonable number of bins for the determination 




(B-v) 



Figure 9. (B — V), V CMD of the cluster with the stars having 
P/j, > 70%. In left panel, we have overplotted the ZAMS given 
by Schmidt-Kaler (1982) while right panel shows the isochrones 
fitting of different ages taken from Girardi et al. (2000). 
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Figure 10. The luminosity function of the stars in NGC 3766 
region with P M higher than 70%. 

of mass function's slope. Fig 1101 shows that the LF contin- 
ues to rise to M v ~ 3 mag and falls precipitously thereafter 
due to incompleteness in the photometry. With complete- 
ness becoming a problem after M v ~ 3 mag (V ~ 16 mag), 
the LF is complete upto V ~ 16 mag. There is no apparent 
structure seen in LF to the level where the photometry is 
complete. 

Using the cluster's parameters derived in this analysis 
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to an increased degree of energy equipartition throughout 
the cluster. The most important result of this process is 
that the higher mass cluster members gradually sink towards 
the cluster center and in the process transfer their kinetic 
energy to the more numerous lower mass stellar component, 
thus leading to mass segregation. The time-scale on which a 
cluster will have lost all traces of its initial condition is well 
represented by its relaxation time Tb. It is given by 



T E 



8.9 x W 5 N 1/2 R s h /2 
< m >V2 log(0AN) 



Log(m) 



where N is the number of cluster members, Rh is the 
half-mass radius of the cluster and < m > is the mean mass 
of the cluster stars (cf. Spitzer & Hart (1971)). The value of 
N and < m > are 283 and 1.5 Mq respectively. Rh has been 
assumed as half of the cluster radius derived by us. Using 
the above relation we estimated the dynamical relaxation 
time Te = 10 Myr for NGC 3766. A comparison of cluster 
age with its relaxation time indicates that the relaxation 
time is smaller than the age of the cluster. Therefore, we 
can conclude that the cluster is dynamically relaxed. 



Figure 11. Mass function for NGC 3766 derived using the stars 
with > 70% and brighter than 16 mag in V. 

and theoretical models given by Girardi et al. (2000) we have 
converted LF to mass function (MF) and the resulting MF is 
shown in Fig llll The MF slope can be derived by using the 
relation log dN/dM = — (l + ^)log(M)+constant, where dN 
represents the number of stars in a mass bin dM with central 
mass M and x is the slope of MF. This relation fit over the 
mass range 11.0 < M/M Q < 1.5 (8.0 < V < 16.0 mag) gives 
a slope of -1.60±0.10 for the MF. Recently, Moitinho et al. 
(1997) derived the MF slope as -1.41±0.08 in the mass range 
12.6 < M/M Q < 2.5. Our derived value of MF slope is in 
agreement with Moitinho et al. (1997). 



8 THE CATALOG 

The catalogue of proper motions with their corresponding 
uncertainties and membership probabilities for ~ 2500 stars 
in the region of NGC 3766 is provided here. In the cata- 
logue, Col. (1) contains the running number; Cols. (2) and 
(3) provide the J2000 equatorial coordinates, while Cols. (4) 
and (5) provide the pixel coordinates x and y. Columns (6) 
to (9) represent relative proper motions and their standard 
errors. Cols. (10) to (13) gives photometric data i.e., B and 
V magnitudes and their corresponding errors. If photome- 
try in a specific band is not available, a flag equal to 99.999 
is set for the magnitude and 0.999 for the error. Cols. (14) 
gives the membership probability P^. A sample of the data 
file is given in Table 



7 MASS SEGREGATION 

To study the cluster dynamical evolution and mass segrega- 
tion, we selected the stars laying within the cluster radius 
and having P M > 70% with V < 16 mag. We subdivided 
the stars into three magnitude range, i.e. 8.0 < V < 11.0, 
11.0 < V < 14.0 and 14.0 < V < 16.0 mag corresponding 
to mass range 11.0 < M/M Q < 7.0, 7.0 < M/M e < 2.0 
and 2.0 < M/Mq < 1.5. In Fig [H] we present cumulative 
radial stellar distribution of stars for different masses. An 
inspection of Fig [12] shows that the cluster exhibit a mass 
segregation effect. Massive stars seem to be centrally con- 
centrated more than the low masses star. To check whether 
these mass distributions represent the same kind of distri- 
bution or not we perform the Kolmogorov-Smirnov (K-S) 
test. This indicates that mass segregation has taken place 
at a confidence level of 97%. Further, it is important to in- 
vestigate that whether existing mass segregation is caused 
by dynamical evolution or the imprint of star formation pro- 
cess. 

One of the possible causes of mass segregation is the 
dynamical evolution of the cluster. Over the lifetime of star 
cluster, encounters between its member stars gradually lead 



9 CONCLUSIONS 

The purpose of this study is to provide a catalogue of precise 
proper motions and membership probability of stars in the 
field of open cluster NGC 3766. We have obtained proper 
motions and astrometric membership probabilities down to 
V ~ 20 mag in 15x30 arcmin 2 area around the open cluster 
NGC 3766. The fundamental parameters, LF, MF and mass 
segregation are derived using the stars of P M > 70%. The 
linear radius of the cluster derived in the present study is 
3.6 pc. The interstellar reddening E(B — V) is found to 
be 0.22±0.05 mag. A distance of 2.5±0.5 kpc and an age 
of ~ 20 Myr have been derived for the cluster. The MF 
slope x = 1.60 ± 0.10 is estimated using the stars brighter 
than 16.0 mag in V. Present analysis indicates that NGC 
3766 is dynamically relaxed and one possible reason of this 
relaxation may be the dynamical evolution of the cluster. 

Finally, we provide the membership probability for 
different types of stars, already found in the literature. 
Here, we also demonstrate that the CCD observations taken 
just 6 years apart can provide accurate proper motions and 
it can be used to separate the cluster members from field 
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Table 3. First few lines of the electronically available catalogue. The stars listed in this table arc found in all the V images. 
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Figure 12. The cumulative radial distribution of stars in various 
mass range. 

stars down to V ~ 20 mag. 
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